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SECTION 1.0 
INTRODUCTION 


One of the secondary objectives of the GEOS- II 
C-Band Systems Project is to study the feasibility of 
using geodetic satellites to both evaluate shipborne 
instrumentation and to determine ship positions in broad 
ocean areas. The purpose of this task is to determine 
whether shipborne C-Band radar tracking, in conjunction 
with ground based tracking, is sufficiently accurate to 
provide instrumentation evaluation and ship position 
estimates. Data from several Apollo tracking ships, in 
particular the USNS Vanguard, was made available for 
this effort. A series of tests, using the USNS Vanguard, 
were carried out in the Port Canaveral and Bahama Acoustic 
Transponder Array Areas. The major portion of the analyses 
reported in this document are the result of preliminary 
investigations using the data from these tests. 

The use of shipborne tracking will be useful in 
several geodetic and geophysical areas provided the reduc- 
tion techniques prove to be feasible and the data of 
geodetic quality. For example, the mapping of areas in 
the vicinity of geodetic anomalies could be performed by 
radar equipped vessels, the location of acoustic trans- 
ponders for use as geodetic reference sources can be 
accomplished in areas not visible from shore, high accu- 
racy navigation updating of position can be accomplished 
for mapping and geophysical experimentation at sea, and 
the use of shipborne trackers as auxiliary or complemen- 
tary stations in geodetic network investigations can be 
accomplished. 
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SECTION 2.0 
TASK OBJECTIVES 


In order to satisfy the overall project objec- 
tives, it has been necessary to establish a set of 
task objectives designed to insure that the information 
flow for the task would be efficient, meaningful and 
useful. Briefly stated the task objectives are: 

a. To understand the data furnished by the 
shipborne radar. This includes the correc- 
tions applied by the ship, the transforma- 
tions performed on the raw data to convert 
it to deck coordinate data, the format of 
the data and what type supporting data is * 
normally furnished by the ship. 

b. To properly calibrate the shipborne radar 
data to remove zero set errors. 

c. To design and implement a set of pre- 
processing procedures to apply the proper 
corrections to the shipborne radar data and 
to provide meaningful navigation system 
information in a "moving tracker" position 
estimation scheme. 

d. To design, checkout and implement a computer 
program capable of estimating the epoch 
position of a moving tracker. 

e. To reduce and analyze actual tracking data 
under a variety of field conditions to 
demonstrate the feasibility of the ship 
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positioning technique and to pinpoint actual 
and/or potential problem areas in the tech- 
nique . 

f. To devise mathematically and physically 
sound techniques for properly accounting 
for sources of error found in e. 

g. To provide a qualitative assessment of the 
technique along with supporting justifica- 
tion for accuracy estimates, and proposed 
further studies. 

The following sections of this report describe 
the approaches used to satisfy these objectives. Problem 
areas are defined in detail and the techniques used to 
overcome these problems are described. The results of 
both simulation and live test analysis and reductions 
are presented with qualitative assessments of the 
results . 
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SECTION 3.0 
METHOD OF APPROACH 


3.1 ANALYSIS OF MOVING TRACKER POSITIONING 

The determination of the position of a ship, or 
other moving tracker, using radar tracking of a satellite 
is very similar in concept to that of positioning land 
based radars. A satellite orbit is determined with the 
help of other stations tracking the same satellite during 
approximately the same time period. Simultaneously, a 
position for the ship is estimated such that the ship 
radar measurements best fit the determined orbit. In 
application, however, there are a number of factors 
which complicate the problem and significantly affect 
the accuracy with which ship position at any time can be 
estimated. 

The following factors distinguish the problem of 
position estimation for a moving radar tracker from 
that for a fixed tracker: 

(1) The variation in ship position from that at 
a reference time must be accounted for at 
some stage of the data processing. 

(2) The difficulty in accounting for short 
period variations in the orientation of the 
radar mount due to ship roll, pitch and 
heading variations makes azimuth and eleva- 
tion angle data essentially useless at 
geodetic satellite ranges for the position 
determinations of interest. Range data is 
thus the only usable radar measurement. 


4 


(3) For most applications, only a limited time 
and number of satellites are available for 
use in the ship position estimation. The 
large number of crossing arcs on all sides 
of the station often available for estimating 
land based stations is conceivably available 
for the positioning of permanent geodetic 
control points. This application is, however, 
somewhat beyond the scope of the present 
report . 

(4) Radar calibration on a ship at sea is more 
difficult to obtain and maintain than for a 
land based radar. Unless very exacting 
calibration and operating procedures are 
followed, range biases of magnitudes varying 
from one satellite pass to the next may thus 
be expected. 

Taken together, all these factors require that the problem 
of precise ship positioning be treated considerably 
differently from land based station positioning. The 
analytical differences will be treated in the following 
two sub-sections. 


3.1.1 Range Measurement Relations for a Moving Tracker 

There are a number of means potentially available 
for determining the position of the VANGUARD at various 
times while at sea. Those for which some data has been 
available on tests to date include: 

1. C-Band radar track of GEOS- II (convertible 
into position) . 
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2. LORAC-B, LORAC-A 

3. SRN-9 (navigation satellite) 

4. Star tracker 

5. SINS, Mark- 19 

6. Acoustic transponder (convertible into position) 

In general, all these methods individually have very 
serious limitations, particularly in accuracy. In 
addition, some are available only at special times, in 
certain geographic locations , or under clear weather 
conditions. Only the SINS system (and, under certain 
conditions, the acoustic transponder positioning 
system) is able to provide a near continuous record of 
ship position. Over a "short" period, its accuracy is 
potentially very accurate. Over a period of many hours, 
its accuracy becomes degraded below that of several 
other methods of positioning. 

The positioning of a ship using C-Band radar data 
alone is an instantaneous fix only if all three measure- 
ment types recorded by the radar (R, A, E) are actually 
used. Biases in the angle measurements, however, make 
this mode of positioning very inaccurate. Such data has had 
very limited utility for results quoted in this report. 

The uses which have been made of angle data have been 
only for transformation of range data and will be briefly 
discussed in connection with the dockside test results 
presented in Section 4.3. 
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Using radar range data without angle data for ship 
position determination required data to be taken at more than 
one time, and required horizontal motion of the ship to be 
accounted for between range measurements. The use of 
the SINS system as the method for obtaining this rela- 
tive ship movement is the logical choice in view of 
its availability, its accuracy, and its continuity. 
Accordingly, we will assume that SINS data can be used 
to relate ship coordinates (<j>, A) at any time t to ship 
coordinates (<})g, Ag) at some earlier time tg . A restric- 
tion must be that t-tg not be too large. For present 
purposes, it will be assumed that t-tg may reach approx- 
imately 2 hours, so that two passes of satellite 
tracking data may be used in the same data reduction, 
with all ship positions accurately related to the ship 
position at the reference time t Q . The desirability of 
using two passes of tracking data is intuitively obvious 
merely from the fact that, in general, one satellite 
track is east of the ship and the other is west of the 
ship. Baring the simultaneous introduction of additional 
error sources, the geometry provided by the two passes 
over and above a single pass cannot but help provide a 
stronger ship position. Nevertheless, single pass 
solutions will not be ruled out of consideration. 

The data reduction for the moving tracker posi- 
tion estimation must then take into account that the 
tracker may be at any position relative to its 
reference time (tg) position, but with the change in 
position available in SINS coordinates which may be 
taken to be latitude, longitude, and height above some 
spheroid A-Aq, H-Hg). In addition, the processing 

must take into account the physical separation of the 
C-band radar from the SINS binnacle to which the SINS 
position is referenced. The latter processing 
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is a transformation of coordinates incorporating both the 
translation and rotation of the true radar coordinates from 
those of the SINS binnacle. The transformation is discussed 
in detail in Appendix A and takes into account the flexure 
of the ship for which measurement data is available. The 
transformation requires the use of SINS heading data but 
makes no use of the radar angle measurements. 

In essence, the SINS and ship flexure data are used to 

relate the position of the radar at any time t to its position 

at the reference time t . Changes in radar coordinates are 

o 

thus not simply differences between SINS coordinates, but 
rather between radar coordinates obtained through trans- 
formations of SINS data. This distinction is of particular 
importance if there is a change in ship heading during a 
track. 

The range measurement made from the moving tracker at 
time t can be expressed functionally as: 

R - ^(X-x) 2 + (Y-y) 2 + (Z-z) 2 (3-1) 

where 

X,Y,Z are ship earth centered fixed coordinates 
at time t. 

x,y,z are satellite earth centered fixed coor- 
dinates at time t. 
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The possible existence of range biases or other such 
systematic errors may be treated in the same way as 
for a fixed tracker and may thus be ignored in the 
present context. In order to make use of the "known" 
relative movement of the ship from its reference 
position, we note that a given (<J>, A, H) is equivalent to 
an (X, Y, Z) with the transformation requiring only 
the spheroid parameters. The relation of the range 
to the initial position of the ship, which is to be 
estimated, is somewhat more involved. Using matrix 
notation, the variation in a computed range measure- 
ment with a change in the reference position of the 
ship may be expressed as 


3R 3X 3 <| » 3(j» n 

6 R = 6X 

3X 34) 3<j> 0 3X Q 


where all symbols on the right hand side now denote 
matrices or vectors. 



3X 3X 3X 


3X q 8X 0 3X q 


3<t> 3 A 3H 


34> 0 3 Aq 3H q 

3X 

3Y 3Y 3Y 

9 *0 

• 

3Y q 9Y 0 3Y q 

3<j> 

3<J> 3 A 3H 

3X 0 

3*0 9A 0 9H 0 


3 Z 3 Z 3 Z 


3Z 0 3Z Q 3Z q 


3<f> 3 A 3H 


J*o 3X 0 3H 0_ 
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I 


3<J> 

3 <)> 

34> 

9< *>0 

3X q 

3H 0 

3X 

3 X 

3X 

^0 

3X q 

3H 0 

3H 

3H 

3H 

3 *0 

3X q 

3H 0 


3R 

3X 


T 



1 

X 

1 

X 

1 


-> 1/R 

Y-y 

, 6X q - 

5Y 0 


Z-z 




(3-3) 


From the assumption that the SINS data provides good 
measurements of the change in latitude, longitude, and 
height, we have the relations 


<J> = <J> 0 + A<|) 

X = X Q + AX ( 3 ' 4 ) 

H = H Q + AH 

With A<J>, AX, and AH constants, the matrix reduces 

to the unit matrix. 
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The other partial derivative matrices are calcu- 
lated from the following relations for transforming 
geodetic coordinates to earth centered fixed coordi- 
nates : 




„ 2 

2 


-1/2 


X = cos 


cos X [H + (1-e 

s in 

<*>) 

] 




„ 2 

2 


-1/2 


Y = cos 

<J> 

sin X [H + a (1-e 
e v 

sin 

<to 

] 

(3-5) 



[H + a e (1-e 2 ) (1-e 2 

2 


-1/2 


Z = sin 

4> 

sin 

<*>) 

] 



where 

a = semi-major axis of the adopted spheroid for 
the earth. 

e = eccentricity of this spheroid. 

9 X 

The matrix yy is obtained by straightforward differen- 
tiation of (2>-5) and will not be given here. The matrix 
9X 0 

Tr-r— is identical in form. 

The computation of the range measurement expressed 
by Eqn. (3-1) and the partial derivative matrices of 
Eqn. (3-2) constitute all necessary calculations that 
are distinct from those which are necessary for esti- 
mating the position of a fixed tracking station using 
range data. In addition, it may be noted that the 
basic mathematical expressions used are identical to 
those for a fixed station. 
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Effects of SINS Errors on Moving Tracker Ranges 


The analysis of the preceding section has been 
based on the assumption that the SINS determination of 
relative latitude, longitude, and height has negligible 
error. Since this assumption is only a convenient, 
and in some cases necessary, approximation, it is highly 
desirable to be able to recognize the presence of SINS 
errors in reduced data. It is also desirable, if at all 
possible, to be able to effectively account for SINS 
errors or to design tests such that SINS errors are 
minimized. 

A brief analysis of the source and form of SINS 
errors is presented in Appendix B for several types of 
ship motion and SINS operation. Some of the conclusions ' 
which may be drawn from this analysis are: 

1. Ship motion during or between radar 
satellite tracks does affect SINS posi- 
tioning errors , but the functional form 
for the errors is not appreciably affected. 

2. Physically distinct error sources produce 
effects which, using data over two 10-15 minute 
segments, separated by a two hour period, 
would be expected to be inseparable on 

the basis of their effects on ship posi- 
tion error. 
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3. SINS position errors, though arising from 
several distinct physical sources, can be 
approximately modeled by the expressions: 


A4> 

II 

> 

h-» 

♦ B.t . 

sin 

(w t + 


AX 

OJ 

< 

It 

* B 2 t ♦ 

C 2 sin 

(wt + 

y 2 ) 


where the A's, B's, C’s, and y's are all 
independent constants. A<f> and AX are 
errors in latitude and longitude, t is 
elapsed time from some arbitrary reference 
time, and w is the Schuler frequency for a 
SINS system operating in an undamped mode. 
When damping is used, the frequency w is 
increased as a function of the degree of 
damping . 


The effects of the SINS error model parameters 
on range measurements enter through the (X,Y,Z) coordi- 
nates of Eqn. (3-1), which in turn depend upon (4>, X, H) . 
Tire SINS error model parameters enter the last set of 
coordinates in a linear manner according to Eqn. (3-6). 
Let Aa be the error in one of the SINS error model 
parameters. Then the range measurement is affected by 
the amount 


3R 3X 3<}) 

6R = 6a (3-7) 

3X 34> 3a 
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The vector ^ and the matrix ^ are given by Eqn. (3-3). 
If we let a be the parameter B^, e.g., the vector 
becomes : 


Partials with respect to the other parameters may be 
obtained using (3-6) in a similar manner. 

For error analyses, which normally require only 
partial derivatives without the need for actual residuals 
the parameter partial derivative expressed by 


3R ax 9<t> 
ax 34) 3a 


is the only quantity needed in order to investigate the 
effects of SINS errors on data reductions which use 
moving tracker data. The manner in which the SINS 
parameter error analysis can be implemented will be 
discussed below. 


3.2 


PROGRAM DEVELOPMENT 


Analysis of moving tracker data used for ship 
position estimation required modifications to two large 
scale computer programs: the A/Omega orbital data 

reduction program and the ORAN orbital error analysis 
program. Although the two programs model the same data 
reduction process, the functions which the programs 
serve are, although complementary, quite distinct. The 
modifications to the two programs are thus completely 
different and will be discussed separately. 


3.2.1 A/Omega Modifications 

The nature of the modifications required for 
A/Omega to incorporate moving tracker data under the 
assumptions made above are of two types: 

1. Input Modifications 

In order for the A/Omega program to be able 
to recognize and use moving tracker data 
properly, the data input format was modified 
to provide: an indicator that the tracker 

was a moving tracker type ; and the changes 
in ship radar position (A<J>, AX, AK) from that at 
some reference time. This reference time 
was basically arbitrary, but was chosen to 
be the time of the first usable data point 
of the satellite pass (or the first satellite 
pass if two successive satellite passes were 
tracked) . 

2. Computations of updated station position. 
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The computation of measurement residuals and 
partial derivatives requires the current position of 
the ship for use in Eqns. (3-1) and (3-3). The A/Omega 
program was originally designed to operate with station 
positions which changed only as a result of a least 
squares adjustment based on satellite tracking data. 
Consequently, the program had to be reconfigured to 
compute a new station position and rotation matrices 
for each data point. Eqns. (3-5) are used to compute 
the new station position and the rotation matrices are 
computed using (3-3). The initial, or reference, 
position of the moving tracker must be retained in core, 
adjusted at the end of each iteration, and used as the 
starting position at the beginning of the next iteration. 

No capability for any SINS error model parameter 
recovery was incorporated into A/Omega, since it was 
not clear that sufficient data existed for the recovery 
of any such parameters. The normal program recoveries 
of measurement biases and station timing errors were 
retained. 


3.2.2 ORAN Modifications 


To investigate the effects of SINS errors on A/Omega 

reductions of moving tracker data, the ORAN program 
was modified for the computation of partial derivatives 

of range measurements with respect to SINS error model 
parameters. The parameters investigated are based on 
Eqn. (3-6) which may be written as: 
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A<f>= (A^ "C^ siny 1 ) +B-^t+ (C^ cosy 1 )sin wt+(-C 1 siny.^) (1-cos ut) 
AX=(A 2 'C 2 siny^)+B 2 t+ (C 2 cosy^)sin wt+(-C 2 siny 2 ) (1-cos wt) 

Since all the parameters are arbitrary anyway, these error 
expressions can be written equally well as 


A<}> = A 1 

♦ B.t + 

C 1 

sin wt 

AX = A 2 

- v - 

C 2 

sin wt 


(1-cos wt) 

(3-10) 

D 2 (1-cos Wt) 


The partial derivatives incorporated into ORAN are the 
partials of range measurements with respect to the 
A, B, C, and D parameters of Eqn. (3-10). The particular 
form of Eqn. (3-10 is chosen such that the constant 
terms will be the position errors at the reference 
time (t=0) . 

The computation of range partial derivatives with 
respect to SINS error model parameters then utilizes 
Eqns . (3-9) and (3-10). Due to the relatively small 

movement of the ship throughout the total tracking 
period, it is possible to make approximations in Eqn. 

(3- 9) which considerably simplify the necessary program 
modifications. In particular, we can set 


3R 3X 3R 3X Q 

3X 3 <j> 3 X q 34> q 


(3-11) 
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Since these partials are already calculated as a part 
of the station position adjustment simulation, the only 
additional computations needed are the type calcu- 

d Ot 

lations based on Eqn. (3-10). These partials are: 


9<f> 

9 A 

3A X 

9A 2 

9<|> 

_____ s 

3A 

3B 1 


9 <J> 

3A 

3C X 

9C^ 

9 <p 

_____ 3 

9 A 

9D X 

9D 2 


1 


t 


sin w t 


1-cos w t 


(3-12 


These expressions are simple functions of elapsed time. 

The ORAN program also required, of course, input 
modifications to request the SINS error model parameters. 
It was then capable of computing the effects of SINS 
errors on: 

1. Adjusted ship range biases. 

2. Adjusted ship position. 

3. Ship range residuals after the least squares 
fit of the range data to the satellite orbit 
adjusting for ship position. 
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The program output was then usable for analyzing the 
A/Omega data reduction results to determine the amount 
and type of SINS errors actually present in the data. 
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SECTION 4.0 
EXPERIMENTAL RESULTS 


The use of shipborne radar data for precise ship 
positioning has been tested in essentially three phases: 

1. Dockside testing to validate the techniques being 
used in the absence of ship motion and with a 
known position of the ship. These tests used 
orbits determined by land based stations and 
tracking geometry similar to that planned for at 
sea tests. The dockside tests thus provided a 
checkout of data handling techniques and gave a 
ship position whose accuracy could be checked. 

The dockside accuracy should be approximately 
equaled at sea except for SINS errors and 
possibly degraded radar calibration. 

2. At sea testing with results checked by analysis of 
ship range residuals and comparison of estimated 
ship position with that obtained from independent 
systems such as LORAC. The results obtained here 
should indicate the accuracy achievable with the 
SINS data for relative positioning accepted with no 
attempt at correction or improvement. Since SINS 
performance can be somewhat variable, some care 
must be taken in extrapolating results based on a 
limited number of tests. 

3. Analysis of at sea tests to determine whether 
effective improvements can be made in the SINS 
data which lead to improved positioning capability. 

For the most part, this analysis consists of comparing 
data reduction results with the ORAN simulations 

of the data reduction, and propagating SINS errors 
into the adjusted parameters and measurement residuals. 
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The results from these three phases should provide 
estimates of ship positioning accuracies currently 
achievable, together with an indication of those pro- 
cedural and data reduction changes most likely to improve 
the results of future tests. 


4.1 DATA SELECTION 

The data selected for the ship positioning 
studies consisted of "raw" range data from one and two 
(successive) passes of GEOS-B. The reasons for these choices 
are given in the following sub-sections. 


4.1.1 Use of Angle Data 

Azimuth and elevation measurements, although made 
by the ship radar, provide position measurements grossly 
inferior to the range measurement. A low angle noise level 
of 10 arc seconds corresponds, for example, to a position 
error of about 100 meters at typical GEOS-B ranges. Biases 
would be expected to be several times this value because of 
the difficulty in angle calibration. In addition, the 
removal of ship motion (pitch and roll) cannot be expected 
to be done without the introduction of appreciable error. On 
the other hand, the radar range noise level is on the order 
of two meters, and its bias should be largely removable 
even for at sea calibrations. Given the exercise of 
sufficient geometry, range data alone is quite capable of 
determining all three coordinates of ship position. Ship 
angle data was therefore not used for ship position studies 
except for some range transformations as indicated in 
Section 4.1.2 and 4.3. 
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4 1.2 Use of Raw Range Data 


In order to be able to provide data for external use 
in fixed tracker data reductions, the ship borne computer 
possesses the capability of transforming (R,A,E) measurements 
to the corresponding measurements for a fixed tracking point. 
This transformation requires two types of input: the location 
of the ship (or radar) relative to the fixed tracking point; 
and the three coordinates of position of the satellite. The 
first of these inputs is provided by the SINS position data. 

The second is provided by the full set of radar measurements 
(R,A,E) smoothed over some interval using filtering appropriate 
for a target in orbit. The transformed data is thus affected 
by SINS errors and radar angle errors (noise and biases) and 
is acceptable data only if there is no better procedure for 
handling it. The use of raw data does not eliminate the effects 
of SINS errors unless there is some auxiliary position measure- 
ment available. It does, however, completely eliminate the 
effects of angle errors since the satellite orbit determined 
by ground based trackers is effectively used as a basis for 
the measurement transformation. For these reasons, ship 
transformed radar data has had very limited use and the 
raw data has been used instead. 


4.1.3 Selection of Data Spans 

In all ship position studies made to date, one or two 
revolutions of satellite tracking data has been used both 
from land based trackers and from the ship. There are a 
number of reasons for both using two revolutions and also 
for limiting the data span to two revolutions. The need for 
more than one revolution of data is because of the much greater 
geometry for the ship which the two passes (rather than one) 
provide, with an associated cancellation of both ship and 
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land based radar error effects. Two revolutions of land 
based data should also provide a somewhat more accurate 
orbit unless a large number of ground stations are 
tracking during the single revolution case (as was never 
the case) . To extend the tracking period beyond two rev- 
olutions means that at least a twelve hour period must be 
used. Orbit errors can be considerably greater for this 
period than for the approximately two hour period needed 
for the two revolutions. In addition, SINS errors be- 
come essentially intolerably large. The one pass results 
are of interest because of the situations in which only one 
pass of data is available, which is the case with SRN-9 
type positioning. 

4.2 ORBIT ERROR EFFECTS ON SHIP POSITIONING 

One of the limitations on the accuracy with which 
ship psoition can be determined using satellite tracking 
data is the accuracy of the satellite orbit itself. With 
the satellite well tracked by highly accurate instruments 
on the ship (e.g., measuring R, A, E) , it would be expected 
that orbit errors would translate one to one into posi- 
tioning errors. Because of the limited amount and type 
of good tracking data, together with the nature of the 
actual orbit errors, there can be a magnification of orbit 
errors when transformed into tracking station error. 

Two station recoveries were simulated using the 
ORAN program for the purpose of estimating the effects of 
various orbit error sources on one and two revolution 
solutions. The orbital solutions were each determined by 
two ground based radars. The one revolution solution used 
data from Wallops and Bermuda and the two revolution solu- 
tion used data from Bermuda and Antigua. The effects of 
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various ground based tracker errors on the recovered 
positions will depend somewhat upon the position of the 
ship relative to each tracking station. The two simula- 
tions performed can, however, be considered as reasonably 
representative of one and two revolution solutions when 
ship range biases must be estimated on each pass. 

Tables 1 and 2 show, the effects of postulated station 
location, range bias, refraction, and geopotential coefficient 
error? on the one and two revolution ship position recoveries. 
For the one pass solution, the largest effects are due to 
ground station biases, for which slightly pessimistic values 
of 3 meters have been assumed. Geopotential coefficient 
errors, as might be expected, have a completely negligible 
effect on this solution. The total uncertainty is approxi- 
mately 12 m, which is comparable to what can be expected 
for the positioning of land based stations using multiple 
passes . 


Table 1 

Effects of Orbit Uncertainty on One 
Revolution (Geos B Rev 7991) Ship Positioning Accuracy 



Effect (Meters) on 

Orbit Error Source 

Latitude 

Longitude 

Station Position \ Lat (5m) 

.1 

1.7 

Uncertainty ( Lon (5m) 

2.0 

6.6 

(Bermuda relative \ 
to Wallops) ) Height (5m) 

1.5 

5.3 

Ship refraction (10%) 

.6 

2.0 

Geopotential Coefficient 
(0.25 APL-SAO) 

.6 

.7 

Bermuda Range Bias (3m) 

1.4 

6.1 

Wallops Range Bias (3m) 

.5 

4.8 

TOTAL 

3.0 

11.5 
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Table 2 


Effects of Orbit Uncertainty on Two 
Revolution (Geos B Revs 8003-8004) Ship Positioning Accuracy 





Effect 

(Meters) On 

Orbit Error Source 


Latitude 

Longitude 

Station Position ) 

Lat 

(5m) 

2.7 


0.2 

Uncertainty 1 

Lon 

(5m) 

2.8 


1.4 

(Bermuda relative 1 

to Antigua) 

* Height (5m) 

6.5 


13.7 

Ship Refraction (101 



2.5 


4.2 

Geopotential Coefficients 
(0.25 APL-SAO) 


2.8 


7.3 

Bermuda Range Bias (3m) 


7.1 


13.0 

Antigua Range Bias (3m) 


.7 


4.3 

RSS Effects 

11.1 

21.2 


The effects of radar biases are still the dominant 
error source for the two pass solution, as shown in Table 2. 
They are actually somewhat greater than for the single pass 
solution, although this may be attributed to less favorable 
geometry for the orbit determination. Geopotential coeffi- 
cient errors now have a non-negligible effect, since such 
errors do produce orbit errors of a few meters for a one 
revolution solution. The total uncertainty of about 24 m 
is considerably greater than that achievably for ground 
based stations with good geometrical coverage. The necessity 
of bias recovery on both passes is largely responsible for 
the magnitude of 24 m. 


25 











4.3 


DOCKSIDE TESTS 


Using two successive GEOS-B tracks by the Vanguard C-band 
radar with the ship in Port Canaveral, it has been possible 
to check out the radar and data handling and reduction 
processes under near ideal conditions. Lack of motion by 
the ship removed the complications of the moving tracker. 

And the position of the tracker was tied to a well surveyed 
point. Position recovery under these conditions should 
give an indication of the accuracy that should be expected 
when the ship motion is perfectly taken in account. 

Table 1 shows the radar position recovered when 
tracks of Revolutions 7972 and 7973 were used in conjunction - 
with land-based tracking by Wallops and Bermuda radars. 

The estimated position, when compared to the survey position 
also shown in Table 3, shows a difference of only 17 meters. 

This figure should be considered a measure of the accuracy 
of the estimated position relative to Wallops and Bermuda. 

For both radar passes, independent radar range 
biases were estimated along with the ship position. Biases 
of -67 and -17 meters were recovered, indicating both 
significant biases and biases which differ between passes. 

The ship radar and land based radar residuals 
(differences between the observed ranges and the calculated 
ranges to the fitted orbit) for Revolutions 7972 and 7973 
are shown in Figures 1 and 2, respectively. The raw radar 
range residuals are denoted by ISTA19. The residuals 
denoted by ISTA66 refer to the low speed data transformed 
by the ship computer as discussed in Section 4.1.2. The 
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trends in these residuals are an indication that something 
was wrong in the transformation process. The ISTA19 
residuals seem quite comparable to the land-based radar 
range residuals. 

In addition to the Vanguard dockside tests 
discussed above, in port reductions have also been made 
of a single revolution of data for an Apollo tracking ship 
in the Sydney, Australia, harbor. The results of this 
test have been reported in Reference 1 and show recovery 
accuracy somewhat less, but still on the same order as 
is indicated by Table 3. 
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TABLE 3 


DOCKSIDE SHIP TEST 
GEOS-B REVS 7972 - 7973 



NAD- 27 POSITIONS 



Adjusted 

Survey 

Latitude 

28° 24’ 31.8” 

28° 24* 31.4 

Longitude 

279° 23’ 44.1" 

279° 23' 44.5 

Height 

7.2 m 

14.6 m 


Orbit Determination by: 

Wallops FPS-16 
Bermuda FPQ-6 
Bermuda FPS-16 
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Figure 1. Residuals for HEOS-R Rov 7Q7? 


4M 



Figure 2. 'Residuals for GEOS-B Rev 7973 


4.4 


AT SEA TEST RESULTS 


At sea Vanguard tracks on three different days have 
been reduced, in all cases using the SINS data for relating 
the position of the ship during track to its position at 
the beginning of track. When tracks on successive revo- 
lutions were available, the initial position used was that 
at the beginning of the first track. 

Results for the three different tests will first be 
given along with some analysis of the recovered ship posi- 
tions and the large systematic measurement residuals which 
were obtained for all at sea tests. In one case, one pass 
of a two pass test will be reduced separately to see how well 
the one and two pass solutions agree. 

It should be noted that a height recovery was not 
attempted for any of the at sea tests. This constraint was 
used because of the limited amount of data on some passes and 
the necessity for recovery of range biases on all passes. 

This necessity was indicated by the dockside tests, which 
showed bias differences on the order of 50 meters from one 
pass to the next (Section 4.3). At sea type data reductions 
were attempted with no bias recovery, with generally disastrous 
results, even when SINS errors were estimated. 


4.4.1 One Pass Solution - Revolution 7991 

Table 4 gives the estimated ship position at the be- 
ginning of track on GEOS-B Revolution 7991. The only avail- 
able comparison position is the SINS position, also listed 
in Table 4. Differences are approximately 440 meters in 
latitude and 30 meters in longitude. A reliable estimate for 
the accuracy of the SINS position is, however, not available. 
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TABLE 4 

AT SEA SHIP TEST 
GEOS-B REV 7991 


Latitude Longitude 

SINS a priori 28° 23' 22.4” 280° 15' 18.0" 

Adjusted (assuming no 

SINS error) 28° 23' 38.8" 280° 15' 24.2" 

Adjusted (including SINS 

error adjustment) 28° 24' 56.0" 280° 16' 35.6" 

Orbit Determination by: 

Wallops FPQ-6 
Wallops FPS-16 
Bermuda FPS-16 
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The residuals for the ship range measurements after 
the ship position (and radar bias) estimation are shown in 
Figure 3. Residuals for the land based radars were some- 
what smaller than the ship residuals and are not shown. 

The SINS velocity error best accounting for the observed 
residuals was estimated using the modified ORAN program. 
Velocity error estimates of 

SV East * 5 ' 5 m/sec 
{V North' 2 ' 8 m/sec 

were found. The effects of this velocity error on the ship 
range residuals are shown as the solid curve of Figure 3. A 
pattern quite close to that of the actual residuals is 
observed . 

The estimated initial position of the ship, taking 
account of the above velocity error, is given in Table 4. 

It will be noted that the ''corrected” position is somewhat 
farther from the a priori SINS position than was the first 
estimated position. In itself, this does not necessarily 
indicate a degradation of the recovered ship position, since 
the accuracy of the SINS initial position is not known. 
However, the estimated velocity error is unrealistically 
large - an indication that the total set of adjusted para- 
meters is somewhat suspect. Further evidence of a weak 
solution is provided by the high correlation (0.99) between 
the two velocity components. The weak solution is due to 
the estimation of 

2 components of ship position 
2 components of ship velocity error 
1 range bias 

using only a limited amount of data. Since the reduction, in 
the same manner, of data from another single pass radar track 
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Figure 3 

Ship Range Residuals for GEOS-B REV 7991 



• Actual residuals 

Effects on residuals 
of estimated SINS 
velocity error 
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produced results which look somewhat better, there is a 
suspicion that the operation of the SINS system during this 
test may have had problems. 


4.4.2 One Pass Solution - Revolution 8004 


GEOS-B was tracked by the Vanguard OBand radar 
and ground stations on Revolutions 8003 and 8004. For 
comparison of results with Revolution 7991, the data for 
8004 was used alone in a ship positioning solution (both 
revolutions of ground tracking were used to determine the 
satellite orbit). The estimated initial position of the 
ship is shown in Table 5. The time of the reference position 
was chosen to be prior to Revolution 8003 so that comparisons 
could be made with the two revolution solution to be discussed 
below. 


In addition to the a priori SINS position for compari- 
son, a LORAC position is also available for this test and 
is listed in Table 5. The agreement of the estimated 
latitude with the LORAC latitude is quite poor, differing 
by about 2500 meters regardless of whether or not SINS 
errors are estimated. Longitude differences are much smaller, 
but are still several hundred meters. The SINS velocity 
error components estimated are 

{V West “ !-2 m/sec 

6V South = 1 '° m/sec - 

As can be seen from the residuals and SINS error effects 
shown in Figure 4, this velocity error provides a reasonably 
good explanation for the systematic trends in the ship range 
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TABLE 5 

AT SEA SHIP TEST 
GEOS- I I REV 80U4 


Latitude Longitude 


SINS a priori 

27° 

7’ 

17.5" 

283° 

39' 

0.6" 

LORAC 

27° 

7' 

32.0" 

283° 

38' 

53.0" 

Adjusted (assuming no 
SINS error) 

27° 

9' 

12.1" 

283° 

39' 

9.2" 

Adjusted (including SINS 
error adjustment) 

27° 

9’ 

7.2" 

283° 

38' 

41.0" 


Orbit Determination by: 

Bermuda FPS-16 
Antigua FPQ-6 
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Ship Range Residuals for GEOS-B REV 8004 
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residuals. The velocity error estimated, however, is 
still too large to be physically acceptable. Again, as 
for Revolution 7991, there are the problems of the adjust- 
ment of five parameters using a limited amount of data and 
a resulting high correlation between the errors in 
the velocity components. 

It should be noted that ship position and velocity 
errors accumulated between revolutions would affect the 
comparisons with LORAC given in Table 5, without really 
affecting the one revolution position determination. 


4.4.3 Two Pass Solution - Revolutions 8003-8004 


GEOS'II was tracked for only a couple of minutes on 
Revolution 8003, far less than the approximately ten minutes 
obtained on some high elevation tracks. Alone, the data 
is of limited utility for ship positioning. In combination 
with a good track of Revolution 8004, the data is quite 
useful because of the additional geometry provided by tracks 
of both sides of the ship. It is necessary, however, for 
SINS errors to be either small or well modeled from the 
beginning of the first revolution track to the end of the 
second revolution track. 

Table 6 shows the adjusted ship position found 
using data from Revolutions 8003 and 8004. Figure 5 
shows the ship range residuals for both revolutions after 
the ship position and bias (one for each pass) estimations. 
Comparison with the LORAC position shows the adjusted 
position to be in error by about 500 m in longitude and 
3000 m in latitude. The range resduals remaining after 
this solution are large (compared to that expected for 
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TABLE 6 

AT SEA SHIP TEST 
GEOS-B REVS 8003 - 8004 


Latitude Longitude 


SINS 

27° 

7' 

17.5" 

283° 

39 ' 

0.6" 

LORAC 

27° 

7' 

32.0" 

283° 

38' 

53.0" 

Adjusted (assuming no 
SINS error) 

27° 

9' 

12.9" 

283° 

39' 

10.9" 

Adjusted (including SINS 
error adjustment) 

27° 

7' 

23.2" 

283° 

39' 

0.6" 


Orbit Determination by: 

Bermuda FPS-16 
Antigua FPQ-6 
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Figure 5 

Ship Range Residuals for GEOS-B REVS 8003 8 8004 
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a two revolution orbital fit for ground stations) and 
trended, particularly for the first revolution. 

The SINS velocity errors providing the best fit 
for the observed residuals are found to be 

6V West ‘ °- 14 m/sec 

SV South = °‘ 53 m/sec ' 

The contributions which these errors make to the actual 
residuals are also shown on Figure 5. Again, the fit appears 
quite good. Removing the effects of these velocity errors 
from the initial ship position, we obtain, as shown in 
Table 6, agreement with the LORAC position to within less 
than 300 meters in latitude and approximately 200 meters 
in longitude. 

The adjusted position, including SINS errors, 
shows closer agreement with the LORAC position than does the 
a priori SINS position. In addition, the SINS velocity 
errors estimated are much closer to physically acceptable 
values than are any of the one revolution solutions. The 
two revolution solution thus appears far superior to the 
one revolution solutions, taking into account the necessity 
for adjusting range biases on each pass and the ship 
initial position and velocity errors. 


The accuracy of the LORAC position was estimated to 
be approximately 18 meters (antenna separation of approxi- 
mately 15.24 meters have not been included in Table 6), 
based on the relative locations of the ship and LORAC trans- 
mitting stations (References 2 and 3) . The system does 
qualify as a good standard in this case. 
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4.4.4 Two Pass Solution - Revolutions 8010-8011 


GEOS-B was tracked on Revolutions 8010 and 8011 by 
the Vanguard radar and two ground based Oband radars. 

The ship passes were somewhat similar to those for 
8003 and 8004, with a short low elevation pass followed 
by a long high elevation pass. Different ground based 
radars were used to obtain the satellite orbit. However, 
the accuracy of the determined orbit should be comparable 
to that for Revolutions 8003 and 8004. The direction 
of the satellite passes (North to South instead of South 
to North) should be insignificant. 

The ship position estimated for the beginning of 
track on Revolution 8010 is given in Table 7 Residuals 
after the adjustment are shown in Figure 6. The residuals 
are somewhat similar to those for the 8003 and 8004 passes, 
with an amplitude approximately twice as large. In addition, 
there is a data dropout during the middle of the second 
pass. Correlations of the data dropout with the data 
quality is unknown. The data reduction assumed the same 
range bias throughout the second pass, but allowed to 
vary from that for the first pass. 

The SINS velocity error best accounting for the 
observed residuals was found to be 


6V 


West 


0.14 


^ North — • 


0 . 


These errors are physically acceptable and do not result 
in very large corrections to the ship initial position 
(compared to the corrections found for Revolution 8003) , 
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TABLE 7 


AT SEA SHIP TEST 
GEOS-B REVS 8010 - 8011 


Latitude 

SINS a priori 27° 6' 41.8" 

Adjusted (assuming 

no SINS errors) 27° 7' 27.1" 

Adjusted (in- 
cluding SINS error 
adjustment) 27° 7' 18.6" 


Longitude 
283° 38' 10.5" 


283° 37' 35.8" 


283° 37' 20.9" 


Orbit Determination by: 


Wallops FPS-16 
Bermuda FPS-16 
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Effects on residuals of 
estimated SINS velocity 


as shown in Table 7. The corrections are less than 

10 seconds in latitude and about 15 seconds in longitude. 

However, the above velocity error does not account very 

well for the actual residuals, as is evident from the 

curves in Figure 6. In particular, the first part of the 

second pass does not fit very well at all, and data 

problems during this period are thus suggested to actually 

exist. No further analysis of the data on this pass 

has been made. Different bias levels on the two parts of 

the second pass is the most probable explanation. However, this 

would require the estimation of one additional parameter 

when the solution is already quite weak. 


4.4.5 Range Bias Estimation 

As has been indicated, range bias recoveries were 
made in all ship test reductions for each. pass. This was 
considered necessary because of the lack of well established 
procedures for precise calibration of the ship radar. The 
biases which were recovered for the various tests are 
listed in Table 8. Values are given for the SINS error 
ignored cases and for the SINS error adjusted cases. 

In absolute terms, some of these biases are quite 
large and inconsistent with even a nominal radar calibration. 
However, the most striking aspect of the Table is the 
large magnitude of the biases for 8003-8004 with SINS 
errors ignored, and the degree to which these values 
are "corrected" when the SINS velocity error adjustment 
is allowed. 
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TABLE 8 


Revolutions 

7991 

8004 

8003 

8004 
8010 
8011 


Estimated Radar Biases 


SINS Error 
Ignored 

SINS Error 
Adjusted 

209 m 

241 m 

-423 m 

-850 m 

2885 m 

135 m 

-404 m 

-91 m- 

-905 m 

-968 m 

192 m 

69 m 
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SECTION 5.0 

CONCLUSIONS AND RECOMMENDATIONS 


The results presented in tire preceding sections 
have shown that C-band radar data from the Vanguard is 
potentially useful for accurate geodetic positioning 
provided the relative ship movement can be accurately" 

taken into account. The results and the implications 
for further analysis may be summarized as follows: 


1. A two satellite pass position estimation 
accuracy on the order of 15 meters is possible 
using a ship C-band radar when relative ship 
position between data points is known from 
other sources. 

2. Ship position estimation using SINS data for 
relative ship motion (Land considered to be 
error free) gives position errors which may 
exceed 3 km. for a two pass solution. 

3. Modeling of SINS errors over a two pass period 
appears possible, with the solution for speed 
and heading errors. 

4. A much stronger station estimation would be 
possible if the radar were well calibrated, 
and no range bias estimation were necessary. 

5. One crucial item reauired in this is approach 
to ship positioning is an independent source of 
relative position. The use of acoustic transpon- 
ders for this purpose should be adequate solution 
under the same conditions. 
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6 . 


A comparison of SINS positions with LORAC 
and acoustic transponder positions should 
help to more clearly identify the magnitude 
and form of SINS errors . 

7 . The preliminary analyses of the Bahama 

experiment should be expanded to incorporate 
the use of acoustic transponder data, and SRN-9 
positions, and future experiments planned util- 
izing the radar calibration information gained 
during this effort. 
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APPENDIX A 

SHIP DATA TRANSFORMATIONS FOR 
MOVING TRACKER IMPLEMENTATION 


Let X be the coordinates of the radar relative to 
o 

a ship origin which should be the center of rotation of 
the ship. This may be closely approximated by the SINS 
binnacle. The coordinate system assumed is X toward 
the bow, Y to starboard, and Z up. 

Let us rotate these coordinates into an (ENV) sys- 
tem which has the X axis east, Y north, and Z up. This 
transformation is for roll, pitch, and heading. 

It is more convenient to derive the inverse trans- 
formation of heading, pitch, and roll. The heading trans- 
formation is a rotation about the Z axis by an amount 
A - it/2 in the negative direction, where A is the heading 
(or azimuth) , 


cos (A - tt/2) 

-sin (A - tt/2) 

0 

sin(A - tt/2) 

cos (A - it/ 2) 

0 

0 

0 

1 _ 

I*” sinA 

cosA 

o“ 

-cosA 

sinA 

0 

0 

0 

1 _ 
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The pitch transformation is then a rotation through an 
angle P about the Y axis in the positive direction 


cosP 


0 


-sinP 


P 


0 


1 


0 


sinP 


0 


cosP 


(A-2) 


Finally, the roll transformation is a rotation about the 
X axis in the positive direction, 


R 



0 



0 

0 


cos rt 

sin 

- sin fi 

cos 


(A-3) 


The transformation from ENV coordinates to local ship 
coordinates is then accomplished by the product of these 
transformations, RPT. Conversely, the inverse transfor- 
mation is accomplished by the transpose of these rota- 
tions, or 


X ENV 


(RPT) T X 

local 


In our case, X local is X Q , 


the radar coordinates. 
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The transformation to ECF -(earth-centered fixed) 
coordinates is made through the relations 


X ECp = (RX) (TX + X ENV ) 

= (RX) [TX + (RPT) T X o 


(A-4) 


These coordinates are, of course, a function of time 
because <J> and X are functions of time. RX and TX are 
given by the equations: 


TX = I -a e 2 sin<|> cos<j> / l-e 2 sin 2 <j> 


a Vl-e 2 sin 2 <J) + H 


(A-5) 


-sinX -sin<(> cosX cos<j> cosX 


RX = 


cosX -sin<j> sinX cos<}> sinX 


(A-6) 


COS(j) 


sin<}) 


where 

a = semi-major axis of the earth 

e 

e - eccentricity of the earth 

<j> = geodetic latitude 

X = geodetic longitude 
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It seems reasonable to assume that the change in 
<J> and X from their values at some epoch time is indepen- 
dent of the initial values of (j> and X. That is, let 


<l> * 4> 0 + A<(> 

X = X + AX (A- 7) 

o 


with A<J> and AX independent of <j> o and X Q . 
Also, we would assume 


H = 


H 


o 


(A- 8) 


since we hardly have any other choice using SINS data. 

For program implementation, however, it is most 
convenient to assume that A<j> and AX are equivalent to a 
change in earth centered fixed coordinates, which is 
independent of <f> Q and X Q . This is the initial assumption 
which should be made. It requires that there be asso- 
ciated with each data point a change in earth centered 
coordinates for the radar. These coordinates may be 
placed on the data cards as follows: 
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Columns 

Format 

Units 

55 - 63 

F9.7 

radians 

64 - 72 

F9.7 

radians 

73 - 80 

F8.3 

meters 



APPENDIX B 

PROPAGATION OF ERRORS IN AN INERTIALLY GUIDED 
SURFACE VEHICLE WITH APPLICATION TO TRACKING 
SHIP ERROR REDUCTION 


SECTION 1.0 

DERIVATION OF SINS ERROR EQUATIONS 


1.1 INTRODUCTION 

In this section, the equations of error propagation 
are developed for a local level inertial navigator. The 
equations are then simplified and the growth of error 
is obtained for the linearized case. The use of vehicle 
trajectory to reduce error growth is presented by com- 
paring the errors obtained for rectilinear versus circu- 
lar motion. 
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1.2 


FUNDAMENTAL ERROR ANALYSIS EQUATIONS 
FOR SURFACE INERTIAL NAVIGATION SYSTEMS 


The following set of equations defines the lin- 
earized dynamics of inertial navigator errors in terms 
of ground-speed velocity and position vectors. The space 
fixed results can be obtained by setting the angular rate 
vectors, ft and p" to zero. 


4» = ^x(ft+p)+e Cl) 

e V = ea - (2ft + p) x eV + eg (.2) 

8 8 

e R = eV - p x eR (3) 

8 

\p - a vector representing the infinitesimal 
rotation from computer to platform axes 

ft - earth rate vector 

p - commanded platform torquing relative to 
earth fixed basis 

- Ground-speed error in earth fixed base 

8 

ea - vector made up of sources of acceleration 
measurement error 
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eg - acceleration error vector caused by error 
in computing g 

t 

R - position vector in lideal platform base 

s - drift rate of platform relative to computer 
axes . 

j 

i 

Equation (1) shows the total drift rate is due to 
applying the total torquring rate (ft+p) about misaligned 
axes plus the contribution of gyro drift. The accelera- 
tion equation and velocity equation are obtained by 
straightforward perturbation of the vector equations 
for ground speed and acceleration. 

Note that the above equation is written in terms 
of the ideal platform analysis which are rotating in 
inertial space at a rate of (ft+p) . 

The acceleration measurement vector can be further 
broken down into the following components. 


ea = b + [K] A t + A T x 

where 

F is the vector of bias terms 

[K] is the diagonal scale factor error matrix 

A^ is measured acceleration. 
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Similarly e can be defined as: 


e = d> + <b . A + 
T o Y iA 


Residual Drift 
<|> 

Y o 


Mass Unbalance 
Drift 

*iA ^SA 


Anisoelastic 

Drift 

^AN 


Scale Factor 
Error 

[C] 


^ SA + ^ AN 


[C] (ft +p ) 


This drift represents the un- 
certainty in the stray gyro 
torques which can not be cali- 
brated out. 

A mass unbalance about the out- 
put axis will produce a torque 
proportional to the accelera- 
tions along the input and spin 
axes. The terms an< ^ *^sA 

refer to drifts caused by mass 
unbalance along the input axis 
and spin axis respectively. 

This gyro drift is produced by 
the unequal deflection of the 
float and is proportional to 
the product of accelerations 
along the spin and input axes. 

The uncertainty in the gyro 
torquing scale factor produces 
this error for gyros that are 
torqued during navigation. [C] 
is assumed diagonal. 
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The most direct method of handling this system of 
linear equations is through the use of the state transi- 
tion matrix concept. The state transition matrix approach 
simplifies the necessary computational bookkeeping and 
presents the equations in a form favorable to accurate 
integration by digital computer methods. A few mathe- 
matical preliminaries are in order: 


X = [A] X (Equations of State) 

[*] = [A] [4> 3 

x (t) = [*(t f t 0 )] X (t o ) 


where equation (4) represents vector shorthand for equa- 
tions (1) through (3) and [<f>] is the state transition ma- 
trix. Note that [A] may be time varying. Equation (6) 
represents the response of the homogeneous equation stated 
in (4) to initial conditions. An appropriate example of 
the state transition matrix approach to inertial naviga- 
tion error analysis is presented in the following formula- 
tion of the equations of state for a local level inertial 
navigator. The equations given utilize the following 
assumptions : 

• Gyro and accelerometer axes are coincident 
with each other and are orthogonal in a 
local level mechanization. 

• Errors caused by instrument misalignment 
and gyro anisoelasticity are negligible. 


• The instrument axes are coincident with 
the ideal local level navigation axes. 

None of these assumptions are particularly restrictive 
and could easily be waived at the expense of a more 
complicated set of equations. 
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The last 18 equations represent the mathematical formalism 
for constant Sources of error. Alternatively, if we view 
the error sources from a stochastic viewpoint, the vari- 
ance of this set of variables is fixed with respect to 
time. 


[ I ] 

nxm 


represents an identity matrix of 
order nxm 


[ 0 ] 
nxm 


represents a null matrix of dimension 
nxm 


[ G ] is the matrix of gravity partial 
differentials . 

The X, Y and Z components given in the matrix 
equation refer to the two level coordinates which are 
arbitrarily oriented and the Z axis which is directed 
along the geodetic zenith (local vertical) . 
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1.3 IMPLEMENTATION OF THE GENERAL SOLUTION FOR 
PROPAGATION OF ERROR 


• Propagation of Error 

The matrix equation given in (5) can be utilized 
directly in the solution for [<j>(t,t o )], where [<j>(t 0 ,t Q )] = 

[I], each element of the <|> matrix representing the re- 
sponse of one state variable to a single initial condi- 
tion. The data required for solution are: 

• Profiles of the nominal trajectory, 

thrust acceleration, velocity and 
position. 

• Estimates (covariance matrix) of the 
initial position, velocity and tilt ' 
errors as well as the gyro and 
accelerometer coefficients. 

The former could be obtained from either a scien- 
tific simulation of the trajectory or from real data 
taken in an operational environment. The latter set of 
variables are dependent upon the uncertainty of cali- 
bration alignment and fix taking as well as statistical 
measurements of instrument performance. 

If these data are available the state transition matrix 
is integrated using the trajectory forcing functions, 
computed from acceleration position and angular rate. 

The estimated values of the error variables are formulated 
as a covariance matrix having diagonal form if the ran- 
dom variables of the system are independent or if the 
interdependence is unknown. The covariance matrix of 
the error is then projected forward in time by the relation: 
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[C ov (X(t),XCt))] = [<*,] [C ov CX(o),XCo)) [<J>] T 

{ 

{ 

In addition, one may use the [<J>] to propagate deter- 
ministic errors forward in time in the form given in 
equation (6). 

• Vertical Channel Stabilization 

As given in the preceeding matrix definitions, the 
navigation error equations are incomplete. In their 
present form, the error equations for vertical are in- 
herently unstable, the vertical channel first order per- 

2 

turbation dynamics having the form of (S - 2g 0 /R Q ) where 
S is the Laplace Operator variable. Therefore, all lo- 
cal level inertial navigation sets have an auxiliary 
sensor and equations to stabilize the vertical velocity 
and position. Normally, this consists of a barometric 
or radiometric altimeter. 

The external data is mixed with the inertial data 
in the form of a linear filter network, the coefficients 
of which are chosen to reduce the error bounds to optimum 
limits. Typically, such an equation would have the form: 


V = a, (R - R ) + f WxWxR dt 
z 1 z z J 


ff A z dt * f G( V 

* a 2 / dt 


dt 
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is the sum of externally derived altitude 
and the earth's radius vector 

compensated value of position component 
in vertical direction 

derivative of R with respect to ideal 
platform basis 

constant coefficients of compensation 
loop 

gravitational acceleration computed 
from external measurements of altitude. 
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SIMPLIFICATION OF THE ERROR EQUATIONS 


If we consider a simplified model of the errors 
due to gravity; 


G x = g o/ R . egx = G x . e 

G = g , , e = G • e R 

y g o/ R g x y y 

G z ' ‘ 2g o/ R ’ eg x = G z .c R 2 


which neglects all terms of order greater than one. Dif 
ferentiating (3) with respect to the ideal frame we have 
after substituting for eV and eV 

o o 


e R + 2W X e R + p X e R + [ (W + fi) • e E] p 
+ (fl 2 - W 2 ) eR + [G] eR = ea. 


If we assume that eR and eR are perfectly controlled 

z z 

by the external device previously discussed and sub- 
stituting W g 2 = g/^ (Schuler frequency) we have for 

the horizontal channel errors: 
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E R x ♦ [(W x t n x ) p x * w 5 2 * si 2 - W 2 ] cR x 


- 2 W 2 ER y ♦ [CW y * Sl y ) P x - P z ] eRy - Ea x 


E Ry * l(Wy ♦ Sly) p y * W s 2 - SI 2 - W 2 ]e Ry 


- 2 We R+ [ (W + ft ) p +p] R = ea 
z x x x' K y K z J x z 


? 

since for all surface vessels, W is between one to two 

s 2 

levels of magnitude greater than W and the 
coupling between the two axes can be seen to be rather 
weak since p and p are quite small, the horizontal 
equation reduces to the following for the Free Azimuth 
(W z = 0) case: 


e R„ + W. 


eR. 


= b + 
x 


K A + 
x x 


K K 


ib A 
T y z 


e R + W, 


eR 


b y * 


K y A y " 


*z A x + + X A z 


In terms of Laplace Transforms, the above equa- 
tions can be rearranged to show the "short term" ( a few 
hours) error response to the most important input types. 


( 8 ) 

( 9 ) 

(S’) 

( 9 ») 
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Response to Gyro Bias Error 


e R (S) = 


zo 


(S 2 + w s 2 ) S 


A (S) 

y 


A • g 
y.O S 


(s 2 + w g 2 ) s 


♦ 


ZO 


A • g 
y xo 6 


e R (S) = 

y 22 

(s z + w s z ) s 


• A fS) + 


CS 2 + w s 2 ) s 


where the cross coupling terms in equation (1) are con- 
sidered negligible and A = [g|. A (S) , A (S) are the 

z x y 

La Place Transforms of the acceleration profiles. 


• Response to Initial Tilt Error 


zo 


e R V (S) = 


S 2 + W 2 

S 




zo 


e R V (S) * 

y ,2 ... 2 


S ♦ W 


• Ry CS) 


y yo 6 


(S 2 + w s 2 ) s 


A • g 
f XO 5 


• A (S) + 


(S 2 + W 2 ) S 


• Response to Accelerometer Bias Error 


S(S 2 + w s 2 ) 


y 22 

SCS^ + W/) 


• Response to Initial Velocity Error 


s 2 4 W s 2 


y 22 

+ w L 

s 
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Gyro Scale Factor Error 




Same form as gyro bias error except that 
4> zo is replaced by C z W z <t> xo is replaced 
by C x W x and <fy Q is replaced by C^W^. 

• Effect of Initial Heading Error 


eR x * V s) ^ZO js 
eR y ’ V x< s > *zo/s 

where 
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1.5 


CHOICE OF TRAJECTORY TO MINIMIZE ERROR 
PROPAGATION 


The simplified equations obtained in the previous 
section show that for constant error sources and step 
acceleration, the error has the form of a sinusoid plus 
a constant, and is therefore, bounded for all errors 
except gyro drift error and gyro scale factor error. 

A brief table of position error for step input error 
sources is given below. 


TABLE I 

ERROR PROPAGATION FOR 
STEP ACCELERATION 


Error Source 



Form of 

Error 

Gyro Bias 

£ 

R *,y 

K io < - t 

i- sin W g t) 

Gyro Scale Factor 

£ 

R *,y 

■ K 20^ - 

sin W g t) 
s * 

Accelerometer Bias 

£ 

R x,y 

- - 

cos W t) 
s 

Initial Tilt 

£ 

R 

*,y 

■ - 

cos W t) 
s 

Velocity 

£ 

V* 

■ - 

cos W t) 
s 

Heading 

£ 

R x,y 

• 

A 

II 

}b • t 

*zo 
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TABLE I (Cont.) 


where 


K 10 ■ «zo • A y>x - 4y,xo • g) / w s ' 


K,„ * (C • W • ftg] + C W • A ) / W ' 
20 v y,x y,y \ ®7 z z y,x' s 


K_ n = b /W ‘ 
30 x,y s 


K 40 ” <J>zo A x,y + ^y,xo • g) / W g 2 


Note that the errors caused by gyro bias and heading 
error are not bounded in time and represent the worst class of 
errors. Consideration of the coupling terms in the ^ equation 
which are ignored here, would worsen the situation 
since in the ¥ x "ip terms the first term is a function 
of velocity over the earth's surface. For the present, 
this effect will be lumped with unforced gyro drift. 

It is clear, in any case, that the application 
of a constant acceleration step is to be avoided since 
the resultant error is unbounded. Therefore, trajectories 
other than straight line motion should be considered if 
vehicle motion is used to minimize system error. 
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Circular Path Motion 


If the vehicle follows a tight circular path at 
fairly low (surface ship) velocity, the acceleration 
profile can be given by 


A (S) 


A y (S) 


AW 


N 


s2 * V 


S 2 * W N 2 


-(A/W n )S 


5 2 * W N 2 


V * 

y 


-A 


s 2 + 


w 


N 


where the ship is assumed pointing at X at t Q , the 

acceleration magnitude is assumed constant and period 

2ir 

of one cycle is . Applying this acceleration pro- 

W N 

file to the equations, we have the interesting re- 
sults of Table II. 


Gyro Bias 


TABLE II 

ERROR PROPAGATION FOR 
SINUSOIDAL ACCELERATION 


eR = K, / — sin W XT t - — sin Wt\ 

11 U N N w s s / 

+ K. ^ (t - — sin W t\ 

12 V w 5 / 
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cos W g t 


Initial Tilt 


tR y - K 11 


W s 2 W N w s 2 


+ cosW xl t 

W N 

W N 


K 12 (* - 7 Sln V) 


eR x “ K 4 iCcosW N t - cosW g t) 


+ K 42 (l - cos W g t) 


W 


N 


eR = K. , (sin W M t sin W t) 

y 41 v N W S 


+ K 42 (1 - cos W g t) 


4>zo A 


Ml 


V - "s' 


M2 


— + f 4> • g) / W 

vr yo,xo BJ / s 


<j> A 
Y zo 


Ml 


2 2 
Wx, - W Z 
N s 


K.- ■ + (<J> • g) /W 2 

42 '• Y yo,xo BJ / s 



Initial Heading Error 


- A *o 

cR - (1 - cos W N t) 

w N 

W N 



sin W XI t . 
N 


A substantial degree of improvement in the error growth 
is shown for the portion of these three terms which are 
contained in the general class of position errors induced 
by faulty heading reference. Clearly, a circular path 
is perferable to a rectilinear one when the location of 
the vessel become secondary to the reduction of error 
growth. Other trajectories could, for a given configura- 
tion, produce more optimum results but more information 
about the inertial navigator would be required before 
the more detailed computation required would be worthwhile. 
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SECTION 2.0 

SHORT PERIOD ERROR ANALYSIS OF A TRACKING SHIP 
INERTIAL NAVIGATOR FOR TWO TYPES OF MANEUVERS 

2.1 INTRODUCTION 

The error analysis equations for the decoupled, 
linearized horizontal channels of a local level inertial 
navigation set are utilized in this section to determine 
the following for two types of between pass maneuvers: 

• Determination of the form of error growth 
for a velocity and position errors in each 
channel 

• Computation of the error propagation produced 
by each error source for a given type of 
maneuver using "typical" magnitudes of error. 

• Comparison of the relative merit of the cruise 
maneuvers . 

The maneuvers referred to in the above describe the 
motion of the tracking ship during a two hour period between 
pass tracking intervals. The two classes of maneuvers con- 
sidered are: 

Case A - Constant Heading, Unaccelerated Motion 

Case B - Circular Motion . 
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ASSUMPTIONS AND APPROXIMATIONS 


The linear uncoupled error equations for the horizontal 
channels of an inertial navigator are utilized throughout. 
The vertical channel error is assumed arbitrarily small 
and uncoupled with respect to the horizontal channels. 

The magnitudes for the error sources used in this report 
are of necessity, somewhat arbitrary, but represent rea- 
sonable estimates of the state-of-the-art instrument error 
budgets and "best guesses" for the initial velocity and 
platform tilt errors. 


Error Source Magnitudes 

Gyro Bias .03°/hr. 

-4 

Accelerometer Bias 1.0 x 10 g 

Initial Tilt 100 sec 

Velocity .2 Kts 

These errors are assumed independent and are the same magni- 
tude for both channels. The effect of gyro torquer error 
and the error induced by torquing the misaligned platform axes 
is lumped in with the residual gyro bias term (.03°/hr.) 
Acceleration induced gyro drifts, accelerometer misalignments 
and all other sources of error are considered second order 
terms for the purpose of this study and are ignored in 
the analysis. * 

The instrument axes (gyro and accelerometer) are assumed to 
be orthogonal collinear, and nominally oriented along the 
computational X-Y axes direction shown below: 


X 

4 


* Y 
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Case A - Unaccelerated Constant Heading 

V NON 10 Kts * 

Heading - Positive Y Direction 


X 

♦ 


* Y 


0 

1 


* 0 

2 


Position(|) represents initial position and is 
the nominal position of the vessel two hours later. 


Error Source 


Gyro-Bias 


Heading Error 


Position Error 


e R„ = 


d> * g 

Y yo 6 


W. 


(t - jL sinW s t) 


eR 


d) • g 
y xo 5 


W 


(t - i- sinW t) 
\ w s / 


eR x " V o ‘ *zo • t 


cR y = 0 
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8 


Initial Tilt Error 


eR x - 


yo 


W 


(i - cos v) 


<p • g 

T xo 6 


eR = 

y 2 

w L 


(l - cos v) 


Accelerometer Bias Error 


eR = 
x 


W 


(i - COSW s tj 


eR y " 2 

W s 


(l • 


Initial Velocity Error 


eV 

x 

eR 3 sinW t 

w s 


eV 

y 

eR = sinW t 

y w s 


For Case A the form of error propagation is then 


R =A»t + B + C sin (W t + <j>) 
x,y s 


where the coefficients A,B,C and depend on the manner in 
which the errors combine together. Note that only 
gyro bias and heading error contribute to the unbounded 
error. A plot of the magnitude of error for the X 
channel is given in Figure 1 for the assumed values of 
drift, bias and tilt given previously. 
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V o 


10 Kts 

• 

<j> 

Y xo,yo 

= 

. 03°/hr 

Y xo,yo 

» 

100 sec 

4> 

Y zo 

• 

^ = 100 sec 

r zo 

b x,y 

= 

1.0 x 10' 4 g 

eV x,y 

s 

• 2 Kts • 


The corresponding equations are given in Table B-l 

TABLE B-l 

X CHANNEL ERRORS FOR CASE A MANEUVERS 
Gyro Bias Error “ 1.83 (t + .41 sin (4.46t)) 

Heading Error = 4.84 x 10 ^ t 

Platform Tilt Error * 1.66 (1 - cos (4.46t)) 

Accelerometer Bias Error * .344 (1 - cos (4.46 t)) 

Initial Velocity Error = .0445 sin (4.46t 


where the error above is expressed in nautical miles 
and t is given in hours . 


B- 26 


PLATFORM TILT ERROR 



TF M Tilt /tfc KIRS 





Case B - Uniform Circular Motion 


NOM 


= 10 Kts 


W. 


N 


" 2ir /. 334 rad/hr 


V NOM COS W N t 


V NOM Sln W N t 


W N V NOM sin V 


W N V NOM cos W N t 


It is apparent that the accelerometer bias and initial 
velocity errors are the same as for the previous case. 
The gyro bias and initial tilt contributions each have 
an additional term produced by the rotary acceleration 
and the heading error term is completely different. 
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Error Source 


Gyro Bias 


Position Error 


-d> • g 

Y X0 6 


w_ 


(t - — sin W t\ 

\ W„ S / 


eR 


W N V NOM ^zo / 1 


V - V 


/ — sinW M t - — sinW t\ 

U N N « s */ 


d> * g 
y xo 6 


W. 


w. 


sinW t 
s 


:R y 


V NOM ^zo 
W N 2 - W s 2 


W N 2 ' W s 2 W N 


\\ 


— cosW t 

w s 

S 


+ — cosW XT t 

w N 

W N 
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The additional terms produced by the circular motion and 
consequent acceleration are given in Table II with numeri- 
cally evaluated coefficients for the given values of error 
terms and speed. A plot of these values over a two hour 
interval is given in Figure B-3. Most notably, the addi- 
tional errors in gyro bias and platform tilt are very 
small in comparison to the gravity induced terms. The 
heading error is reduced by an order of magnitude for the 
circular motion case, but for the velocity and heading 
error magnitudes considered, this term is quite small. 


Thus, unless the heading error magnitude is much greater 
than that utilized here. The growth of error is essentially 
the same for the two types of maneuvers considered. 

If, however, a significant degree of unanticipated ship’s 
acceleration occurs the circular maneuver would be superior 
since some cancelation of the effects of the pulses can 
be expected. For each acceleration impulse of magnitude 
A, the resultant additional error terms are: 


Gyro Bias 

cR" 

x,y 

-e-* 

N 

o 

> 

'C 

X 

|l - cosW g t 

w 2 

s 

Platform Tilt 

eR" 

x,y 

<|> A 

zo y,x 

sinW t 
s 

W 

s 
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Platform Tilt 





cosWgtj 


eR 


x 


+ 


+ <j> W KT V xinM 
zo N NOM 

"n 2 - W s 2 


cosW^t - cosW g t 


+xo • * 


w. 


1 - cosWgt 


eR 


•d> • g 

y xo 6 

W N 2 - W s 2 


/ WN \ 

( sin V ‘ — sinW s t ) 


Heading Error 


'NOM 


eR 


W 


N 


ill /l - cosW 
r zo \ 




'NOM 


eR 


y W 


N 


+zo sinW N t * 


B-31 



TABLE B-2 


CIRCULAR MOTION INDUCED ERRORS 
IN X AND Y CHANNELS 


Gyro Bias Error 

eR' 

X 

= 1.58 xlO ^ sin(18.8t) 

- 6.62 x 10" 5 sin (4.46t) 


eR' 

= 1.49 xlO" 5 - 6.68 xlO -5 


y 

cos (4 . 46t) 



+ .84 xlO" 5 cos (18 . 8t) 

Platform Tilt Error 

eR' 

* 2.73 x 10' 4 (cos (18 . 8t) - 


X 

cos (4 . 46t) ) 


eR' 

= 2.73 x 10 -4 (sin(18.8t) - 


y 

4.42 sin(4.46t)) 

Heading Error 


= 2.58 xlO' 4 (1 - cos (18 . 8t) 


£R y 

= 2.58 xlO -4 sin (18. 8t) 


For eR 


in Nautical Miles and t in Hours . 



FLAT FORM TILT l» 




SECTION 3.0 

SHORT PERIOD ERRORS FOR A 
DAMPED INERTIAL NAVIGATOR 


3.1 INTRODUCTION 


The damping of a marine inertial navigation system by means 
of an external source is generally accomplished by subtracting 
the reference velocity from the SINS velocity and adding the 
filtered difference to the measured acceleration as a correction. 
In terms of the short period error model, this amounts to provid- 
ing negative feedback, i.e., damping for the error term at the 
expense of introducing a new error source. The filter used for 
the damping is considered a gain in the following which is con- 
sistent with the level of approximation implicit in the short 
term error model. 


3.2 DAMPED INERTIAL NAVIGATOR ERROR EQUATIONS 


« *x * V R x * W s £R x 


b x + K x A x 


+ \p A - \p A 
r z y r y z 


+ K, <5V 
1 xr 


cR y + K l£ R y + W^eR y 


b + K A 

y y y 


*z A y + Vz 


+ K, <$V 
1 yr 
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The short period errors produced by the damped system 
are shown in the following for the x channel in Figure B-4 
and in the equations below based upon the following assump- 
tions and approximations. 1 

1) The vehicle is a slowly moving local level system, 
such that the horizontal (x,y) channels are 
uncoupled and independent. 

2) For the time interval of interest the effect of 
the earth's rotation is negligible. 



i|> v - platform tilt about y axis e R - Position error 

e y ‘ gyro drift about y axis . - Damping gain 

6V^ X - Velocity reference error 
in x direction 


FIGURE B-4; X CHANNEL ERROR MODEL 
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X Channel Errors 


• GYRO Drift 


t _ h + h e '(% n n 


2 2 

w: wf 

s s 


cos (Wgt) 


• Platform Tilt 


• Velocity Error 


(1 - — ) / 


2 , " (It l/2)t 

— -) / w sin W t 


CR X = R 4- y 


-(K 1/7 )t/ K 

1 - e 1/2 ( cos (W't) + — sin (W" t) 

I s 2W . s 


K 1 6V xr 


|l - e ( 1/2 ' ) ^cos (W^t) + -i- sin (W^ t) jj 
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• Accelerometer Bias 



b x 

/ K l/2^ t , K 1 

\ 1 - e 1 cos (ICt) + 


W s 

\ ' s 2Wg 

% 

where 

W s 

g/R Schuler frequency 


W' 

S 

K 2 

(W 2 - — ) 1/2 Reduced Schuler 
s 4 

- 

e 

y 

Y gyro drift rate 


*y 

Platform tilt about Y axes 


« v xr 

Velocity reference error 


b x 

Accelerometer bias error 


t 

time in hours . 


sin (W't) 


frequency 
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I 


In Marine Systems the design coefficient of damping is 
typically between .3 and .7, and since: 


K 1 = 2 C W s 


Let X, = .4 and K, = . 81V 

Is 


W =4.46 (in hours units) . 
s 

3.3 RESIDUAL ERRORS 

The damped inertial system will, .for step errors, be 
well damped (to 5 percent) in three time constants, i.e., 

=sl 2/3 hrs. assuming that no other errors are introduced the 
residual errors are the steady state and/or ramp errors associated 
with the damped SINS and represent the ultimate limit of accuracy 
for a system of invariant error sources with no external aids, 
i.e., without any fix information. 


• GYRO Bias 

• Platform Tilt 

• Velocity Reference Error 

• Accelerometer Bias 


.06 (+ .177) 


naut. mi. 
.001' o/hr. 


1 naut. mi./-rr- 
min 

.177 naut. mi. /knot 
- 4 

.3 naut. mi./lO g 


« 
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I 


These values are useful in determining the upper bounds 
of the error magnitudes permissible during the tracking period. 
For example, suppose the mission requires 100 meter accuracy 
during a two-hour tracking period and given that the system 
is perfectly quiescent at the initiation of tracking, we can 
determine the single source error bounds as follows: 


Resultant Position 
Error After Two Hours 


Permissible Single Error 
Source Magnitude 


GYRO Drift of .000385 °/hr. 


or 


Velocity Reference Error of .28 KNOTS 


or 


Platform tilt of 3 sec 


or 


^ Accelerometer Bias of 1.67 x 10"^ g • 

In terms of error sources number two and three we are 
proscribing step changes in EM Log error or drift current of 
greater than .3 knots for the former and deflection of the 
local vertical step changes of more than 3 seconds of arc for 
the latter during the two-hour period. Since we have considered 
only one of two channels, the allowable error is even less 
than shown, however the magnitudes given are good ’’ball park 
numbers" and serve to illustrate that according to SINS 
specifications the ability of the unaided SINS to maintain 
100 meters accuracy over a two-hour period is -at best marginal. 


100 meters 
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3.4 AIDED INERTIAL OPERATION (SKOR) 


The Sperry Kalman Optimal Reset (SKOR) program can be 
utilized during the tracking period in the STORE Mode to 
compute the best estimates of system errors obtained by 
performing a Kalman filtering operation utilizing the Marine 
Star Tracker (MST) measurements. 

The output of the SKOR program is printed out at six 
minute intervals and should provide a considerable improvement 
over the unaided inertial system. However, the required usaqe 
of the MST is not possible if there is much cloud cover, 
in which case the system accuracy reverts back to that of the 
unaided SINS discussed previously. 



